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Since the beginning of this century we have attended a blooming of the gravity-mode research thanks to the unprecedented 
quality of the data available, either from space with SoHO, or from the ground-based networks as BiSON or GONG. From 
the first upper limit of the gravity-mode amplitudes fixed at 10 mm/s at 200 /iHz given by Appourchaux et al. (2000), 
on one hand, a peak was supposed to be a component of the £=\, n=\ mixed mode (Garcia et al. 2001a, b; Gabriel et 
al. 2002) and, on the other hand, a couple of patterns -multiplets- were attributed to gravity modes (Turck-Chieze et al. 
2004; Mathur et al. 2007). One of these patterns, found around 220 /iHz, could be labeled as the £=2, n =-3 g mode, which 
is expected to be the one with the highest surface amplitude (Cox and Guzik 2004). Finally, in 2007, Garcia et al. were 
able to measure the fingertips of the dipole gravity modes looking for their asymptotic properties. In the present paper 
we present an update of the recent developments on this subject with special attention to the 220 /iHz region, the dipole 
asymptotic properties and the impact of the incoming g-mode observations on the knowledge of the solar structure and 
rotation profile. 
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1 Introduction 

In the early 80's, right after the beginning of the helioseis- 
mology, the gravity-mode (g modes) searches started. Sev- 
eral groups looked for both individual modes and the sig- 
nature of their asymptotic properties (see for example: De- 
lache and Scherrer (1983); Frohlich and Delache (1984); 
Isaak et al. (1984); Palle and Roca-Cortes (1988) and the 
reviews by Hill et al. (1991) and Palle (1991)). Unfortu- 
nately, none of these candidates could be confirmed as grav- 
ity modes by more recent observations. In the middle of 
the 90's, two new observational facilities became fully op- 
erational in helioseismology: the complete deployment of 
the ground-based networks (e.g. BISON 1 and GONG 2 ) and 
the launch of SoHO 3 (with three helioseismic experiments: 
GOLF 4 , VIRGO 5 and SOI/MDI 6 ) expanded the research 
opportunities. 

With the passage to the new millennium, we attended 
a blooming of g-mode research based on the quality and 
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1 Birmingham Solar Oscillation network (Chaplin et al. 1996) 

2 Global Oscillation Network Group (Harvey et al. 1996) 

3 Solar and Heliospheric Observatory (Domingo et al. 1995) 

4 Global Oscillations at Low Frequency (Gabriel et al. 1995) 

5 Variability of solar IRradiance and Gravity Oscillations (Frolich et al. 
1995) 

6 Solar Oscillations Investigation/Michelson Doppler Imager (Scherrer 
et al. 1995) 



accumulation of data. The first efforts were concentrated 
in looking for individual modes. Indeed, in 2000, Appour- 
chaux et al. looked for individual spikes in the power spec- 
trum of several instruments (VIRGO, BiSON and GONG) 
concluding that, with 90% confidence level, no signals were 
detected. Although they could not identify any g-mode sig- 
nature, an upper limit of their amplitudes was established: at 
200 /iHz, they would be below 10 mms -1 in velocity, and 
below 0.5 parts per million in intensity. The same year, us- 
ing the GOLF instrument and different power spectrum es- 
timators, Garcia et al. (2001a, 2001b) found a peak around 
285 /iHz that could be interpreted as one component of the 
£ = l,n = 1 mixed mode (up to 98% confidence level). 
Later, in 2002, Gabriel et al. (2002), using a similar statis- 
tical approach, confirmed the existence of this mixed-mode 
candidate with a longer data set of the same instrument. Be- 
sides, two other structures at around 220 and 252 /iHz - 
already studied with 2 years of GOLF data and reviewed 
by Gabriel et al. (1999) - were highlighted to be potentially 
interesting due to their persistency with time. 

In order to reduce the threshold while maintaining the 
same confidence level, Turck-Chieze et al. (2004) looked 
for multiplets instead of spikes and applied it to the g-mode 
research with GOLF. This study led to several patterns that 
were considered as g-mode candidates. Unfortunately, this 
analysis fails to properly label these modes with the cor- 
rect degree and orders (n, £, m). In particular, the struc- 
ture around 220 /iHz retained our attention. In fact, Cox and 
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Guzik (2004) deduced theoretically that the g mode with the 
highest surface amplitude would be the £ = 2, n = —3 ex- 
pected at 222.14 /iHz (from the Saclay seismic model com- 
puted by Couvidat, Turck-Chieze and Kosovichev (2003)). 
Although the precise structure of this pattern changed with 
time, the main features still remained above the statistical 
threshold when 3000 days were used (Mathur et al. 2007 
and references therein). 

In 2007(a), instead of looking for individual g modes, 
Garcia et al. found the nearly constant spacing in period of 
the dipole g modes between two consecutive modes of the 
same degree (£=l) and consecutive radial order n using 9.5 
years of GOLF data. 

Today, we face a situation where the fingertips of the 
gravity modes have been found collectively as well as some 
individual g-mode candidates. However it would be extre- 
mely important if both discoveries could be confirmed by 
other instruments and also with longer data sets. Thus, in 
section 2 we extend the GOLF data series, on one hand, by 
more than 40% since the last time we have studied the g- 
mode candidate around 220 /iHz and, on the other hand, by 
more than 20% for the study of the asymptotic properties 
of the dipole modes. Then, in section 3, we study -using 
different instruments- the region of the g-mode candidate 
£ = 2, n = — 3 and the behavior with time of one of the 
multiplet components at 220.7 /iHz. We conclude in section 
4 by doing some theoretical prospects on how the studies of 
the physics and dynamics of the solar interior could bene- 
fit in the near future by adding physical constraints coming 
from the detection of g modes. 

2 GOLF 11.5 year-long series 

Nowadays, we have cumulated nearly 11.5 years (4182 d) 
of GOLF data from 1995 April 11 to 2007 September 22. 
These data correspond to time series calibrated into velocity 
following Garcia et al. (2005). As we have seen in the pre- 
vious section, two different strategies have been followed 
to analyze the data: to look for individual modes above 150 
/iHz, or to look for the asymptotic properties between 25 
and 140 /iHz. 

2.1 Looking for individual multiplets: 220 /iHz region 

The most interesting case found in the previous analysis of 
the GOLF data -when individual multiplets were searched- 
was the region around 220 /iHz. Indeed, an interesting pat- 
tern of peaks has been tracked with time since the very be- 
ginning of the mission when only two year-long data sets 
were available (Garcia and Turck-Chieze 1997; Palle and 
Garcia 1997). These analyses were based on the recurrence 
with time of some peak-patterns and their findings were 
summarized in Gabriel et al. (1999). 

Later, this region was studied with longer datasets. In 
particular with 1290 and 2975 days long (see Mathur et 
al. 2007 and references there in). First, with 1290 days and 



using the average of 4 tapered periodograms, a triplet was 
found centered at 220.72 /iHz and with two components at 
each side at: 220.10 and 221.28 /iHz. Another peak raised 
above the 98% confidence level at 218.95 /iHz which, at 
that time, seemed quite far to belong to the same mode. 
However, all together could be detected as a quadruplet with 
more than 98% confidence level in a 10 /iHz window. Later, 
with 2975 days and using the average of 8 tapered peri- 
odograms, the same peaks stand above the statistical thresh- 
old and a new one (at 219.59 /iHz) raised between the left 
peak and the others forming a quintuplet as a whole. 
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Fig. 1 Power spectrum density of 4182 day-long GOLF veloc- 
ity time series. The vertical dotted lines are the components of the 
multiplets previously seen with GOLF. The vertical dashed lines 
are the expected theoretical frequencies for the g modes £=2 and 3 
in this frequency region from the Saclay seismic model. The hori- 
zontal line is the 90% probability threshold for a quadruplet. 

Today, the same type of periodogram estimator has been 
computed using 4182 day-long series (see Figure 1). To 
keep the same significance in the statistical analysis we have 
maintained the same number of points in the studied win- 
dows. Therefore, due to the increase in the power spectrum 
resolution, we were obliged to reduce the frequency win- 
dow from 10 to 7.1 /iHz. The resultant spectrum is closer 
to the situation we had with 1290 days compared to the one 
we had with 2975 days series. In fact, the peak that raised in 
between the two groups (at 219.59 /iHz) is now embedded 
into noise and the one at 220. 10 /iHz is at the detection limit. 
Therefore a quadruplet is still detectable if we consider the 
uncertainty related to the computation of the 90% threshold 
(see the discussion on Turck-Chieze et al. 2004). The theo- 
retical frequencies for the g modes £ =2, n = -3 and £ = 3, 
n = -5, calculated with the seismic model are superimposed 
on Figure 1 . A study on the uncertainty on these predictions 
has shown that depending on the physics included in the so- 
lar models (abundances, microscopic diffusion, horizontal 
diffusion in the tachocline), we can have an uncertainty on 
the frequencies of 3-4 /iHz in this frequency range (Mathur 
et al. 2007). 
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Fig. 2 Power Spectrum of the GOLF power spectrum density computed between 25 and 140 /xHz. Top panel: as a function of the 
period; bottom panel: as a function of frequency. The shaded regions are the APi and AP2, m = regions and their first two harmonics 
(blue and orange respectively). 



2.2 The dipole asymptotic properties 

Using 3481 days of GOLF time series, Garcia et al. (2007) 
uncovered the presence of a pattern in the power spectrum 
(RS.) of the power spectrum density between 25 and 140 
/iHz that could be interpreted as the quasi constant separa- 
tion in period between consecutive radial orders n of dipole 
(1=1) modes. This separation (APi) was characterized as 
a pattern of peaks with a maximum of 6.5 a and a total 
power of 2.9 times more than the average power of the rest 
of the power spectrum, a is the standard deviation defined 
as (1/(7V - 1) EjLo^i - x) 2 ) 1 / 2 , N is the number of 
points between 4 and 50 min in the periodogram, and Xj are 
the individual bins. 

With the new series 20% longer (see Figure 2), the struc- 
ture has grown up to a maximum power of 8.6 a and a total 
power of 4. 1 times higher than the average power of the P.S. 
It is also important to note that the highest peak of the spec- 
trum between 4 and 50 minutes is the one attributed to the 
APi structure. Besides, the other highest peaks correspond 
to the regions where we expect power from the harmonics 
of this peak as well as the ones corresponding to the period- 
icity, AP 2 of the quadrupole modes (£=2). In Figure 2, we 



have only shaded the regions where the AP 2 for the m=0 
components and its harmonics are expected: they are the 
only m components, the zonal ones, that are not modified 
by the rotational splitting. Indeed we see near the region of 
the second harmonic two high peaks at ~ 0.003 15 Hz and ~ 
0.00385 Hz that could be the signature of the m — ±2 com- 
ponents. A deeper study will be done in the future to explore 
the structure of the AP 2 periodicity and its harmonics. It is 
important to remember that the position of the shaded re- 
gions are almost insensitive to the choice of the solar model 
(Mathur et al. 2007). However, their sizes depend rather on 
the simulated rotation rate, mainly for the AP 2 and this is 
the reason why we have only considered the zonal compo- 
nents. In the case of the APi the resolution of the P.S. is not 
enough to correctly separate the components even for very 
high core rotation rates. The two-peak structure observed 
in the region between 22 and 26 minutes could be due to 
the high level of noise rather than the signature of the m ± 1 
components. Further studies would also be necessary before 
any claim in this direction. 

The conclusions obtained after the analysis of the recon- 
structed waves built after filtering the APi region between 



www.an-journal.org 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



4 



Garcia et al.: Update on g-mode research 



22 and 26 minutes and its first harmonic are qualitatively the 
same as in the previous shorter study (Garcia et al. 2007). 



3 The region around 220 /iHz. 

In the previous section we have seen that there is a statis- 
tically significant (above 90% confidence level) structure 
around 220 /iHz using GOLF data and this structure -which 
evolves with time- has been present in these data continu- 
ously during the last 11.5 years. In adition, when using the 
longest data set available for the SPM of the VIRGO exper- 
iment (4098 day), this feature also appears and with higher 
significance than in GOLF. So, in this section both data sets 
will be used to analyze the time evolution details of this fea- 
ture. Figure 3 shows the power spectrum density in the same 
frequency region than in the GOLF case (see Figure 1). It 
corresponds to 4098 day-long time series of the blue SPM 
channel from 1996 April 11 until 2007 June 30. Compared 
to the multiplet components detected with GOLF, only the 
one at 220.7 /iHz component is clearly visible with a high 
signal-to-noise ratio (> 8 a). Some power excess is also 
present at the most right peak at 221.28 /iHz but slightly 
shifted compared with the frequency extracted from GOLF. 

Besides those components, some other peaks are present 
in both spectra. In particular, around 218.2 and 217 /iHz, 
there are some patterns in common in both instruments with 
peaks in a range of 6 to 8 a. Therefore, this frequency re- 
gion is of particular interest, deserving the use of all the best 
available data sets. 
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Fig. 3 Power spectrum density of 4098 days of intensity time 
series computed from the SPM blue channel. The vertical dotted 
lines are the components of the multiplets previously seen with 
GOLF. The vertical dashed lines are the expected theoretical fre- 
quencies for the g modes £=2 and 3 in this frequency region from 
the Saclay seismic model. The horizontal line is the 90% probabil- 
ity threshold for a quadruplet. 



3.1 Time evolution as seen by VIRGO instruments 

The VIRGO experiment and its performance has been al- 
ready reported in detail (Frohlich et al. 1995; Frohlich et 
al. 1997). One of its instruments is a three channel full-Sun 
PhotoMeter (SPM) which measures solar irradiance through 
5 nm wide filters centred at 402 nm (blue), 500 nm (green), 
and 862 nm (red) that look at the Sun as a star with a 60 s ca- 
dence. VIRGO experiment also contains two different acti- 
ve-cavity radiometers (DIARAD and PMO) for monitoring 
the solar "constant" and a Luminosity Oscillation Imager 
(LOI) with low resolution (12 pixels) for the measurements 
of the radiance distribution over the solar disk at 500 nm. 

From a VIRGO/SPM time series of 3620 days (from 
1996 April 1 1 to 2006 March 3), 57 time series of 800 days, 
each one shifted by 50 days consecutively, have been ex- 
tracted and the respective power spectra computed. In this 
way, the temporal evolution of the spectral signals can be 
studied. The power spectra of the 800 days time series are 
computed with a zero padding factor of 5, thus increasing 
the number of bins by the same factor. The resultant fre- 
quency resolution is 0.00289 /iHz but the points are corre- 
lated and not statistically independent (see the discussion in 
Gabriel et al. (2002) about the use of zero padded spectra). 

In the VIRGO/SPM data a clear and stable signal at 
220.7 /iHz is observed during the 10 years of the SoHO 
mission. The signal is well observed in all three channels of 
SPM (blue, green and red) as shown in the Figures 5 (1, 2 
and 3). The corresponding average power spectrum of the 
57 power spectra has also been computed for each chan- 
nel and a sample of them is shown in Figure 4. In order 
to look if this signal is also observed in the others VIRGO 
instruments, the same analysis have been undertaken using 
the other data instruments: LOI imager (averaging the 12 
pixels) and the PMO and DIARAD active radiometers, pro- 
viding positive results on the presence of the same feature 
in all of them (see Figures 5 (4, 5 and 6)). 
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Fig. 4 Average of 57 power spectra corresponding to the blue 
channel of VIRGO/SPM (figure 5) where the signal at 220.7 /iHz 
is clearly visible. 
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The possibility that this signal could be induced by some 
instrumental effects has also been investigated. In particular, 
temperature fluctuations of the detector. Each of the three 
SPM channels measurements are corrected by a quantity 
proportional to their respective value of the detector tem- 
perature. This correction is applied at the level 1 stage of 
the data pipeline, so data used in the present analysis are 
already multiplied for this quantity. This correction is: 

SPM C 

riannei — (1 + 0= riannei 

(TS C 

riannei 

- 293.15)), 

where channel means blue, green or red; C channel is a con- 
stant for each channel and TS 'channel is the temperature of 
each of the three detectors. 

In Figure 6 we show the SPM blue signal (1), the SPM 
blue detector temperature fluctuations (2) and in (3) a zoom 
of the fainter region of (2) (time series 20 to 57). We do not 
see any correlation between signal and temperature being 
the temperature fluctuations almost three order of magni- 
tude smaller than the 220.7 /iHz signal. 

3.2 Time evolution as seen by velocity instruments 



4.1 Constraints on the rotation profile 

The rotation of the solar interior is well known down to the 
frontier the solar core (e.g. see Thompson et al. 2003 and 
references therein). Unfortunately, p modes do not contain 
enough information to properly constrain the rotation inside 
the nuclear core (e.g. see Couvidat et al. 2003; Garcia et 
al. 2004), even using high radial order, low-degree p modes 
(Garcia et al., 2008). The analysis of the APi structure and 
the comparison with solar models (with different artificial 
solar rotation profiles inside the core) favored -in average- 
a faster rotation rate in the deepest layers. However it was 
impossible to constrain with accuracy the profile. 

To progress in our knowledge of the inner core we need 
to include information coming from the analysis of low- 
degree low-order p modes, mixed and individual g modes 
(e.g. see Provost, Berthomieu and Morel, 2000). Recently, 
Mathur et al. (2008) have shown the influence on the in- 
ferred rotation rate of the solar core when 1 and 8 simu- 
lated g modes are introduced in the 2D inversions based on 
the Regularized Least-Squares method (for more details see 
Eff-Darwich and Perez Hernandez, 1997). 



This signal at 220.7 /iHz has a high visibility in the inten- 
sity data of the VIRGO experiment. It is also interesting to 
verify the behavior of this region in velocity. A preliminary 
analysis of GOLF, MDI (integrated data) and GONG 7 (sun- 
as-a-star data) -ground based Network- has also been done 
in the same way as for VIRGO. Figure 6 shows the resultd 
for GOLF (4), MDI (5) and GONG (6). The visibility of this 
signal in velocity seems to be much lower than in intensity 
but a trace of the signal is observed mainly in GOLF and in 
GONG data. 



4 Inferring properties of the radiative 
interior 

During the last years we have seen the detection of some 
patterns that have been considered as g-mode candidates. In 
this paper we have seen that several instruments have inter- 
esting features in the power spectrum. Some of them have 
been presented in the spectrum for the last 4000 days and 
the careful analyses of the VIRGO-SPM house-keeping pa- 
rameters did not show any similar instrumental behavior. 
Everything seems to point to a solar origin. Whether it is 
signal from the g modes or just a very unlikely combination 
of granulation noise, the future will tell us. The possibility 
of measuring individual g modes could be just in front of us. 
It seems that we have likely measured power coming from 
the cumulative effect of several of such modes. Therefore, 
we should be prepared to incorporate all this information 
to better constrain the structure and the dynamics inside the 
radiative zone and, in particular, in the core. 



7 Global Oscillation Network Group (Harvey et al., 1996) 
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Fig. 7 Equatorial rotation rate below 0.4 Rq inferred from the 
use of p modes only, when 1 g mode is introduced and when 8 g 
modes are introduced. The simulated rotation profile is the contin- 
uous blue line. 

Figure 7 summarizes the conclusions of this work. In- 
deed, an artificial profile has been computed in the core: 
a step profile having a rate of 350 nHz in the 0.1-0.2 Rq 
region and a rate 3 times larger than the rest of the radia- 
tive zone below 0.1 Rq. Although this profile is unlikely to 
be realistic, it enables us to characterize the quality of the 
inversions, as the code has difficulty in reproducing these 
steep gradients. Using this profile, the rotational splittings 
of the p and g modes have been computed and the inver- 
sion calculated with three different sets of modes. The first 
set contains only the same p modes than those available in 
the real case. In fact, the error bars have been taken from 
real observations. This is the reference inversion (continu- 
ous line with errors in Figure 7). However this result is reli- 
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Fig. 5 Time evolution of the 220.7/iHz signal as observed in all instruments of the VIRGO package: SPM blue, green and red; LOI, 
PMO and DIARAD instruments. Each vertical column corresponds to the PSD of 800 days shifted by 50 days. 
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Fig. 6 Time evolution of the 220.7/iHz signal for SPM/blue (1); detector temperature fluctuations for SPM/blue (2); (3) is a zoom of 
(2). In (4), (5) and (6) the same analysis but for velocity data from GOLF, MDI and GONG instruments. Although weaker in velocity 
than in intensity, a trace is observed mainly in GOLF and GONG data. Each vertical column corresponds to the PSD of 800 days shifted 
by 50 days. 
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able down to 0. 16 or 0.2 Rq depending on the set of p mode 
used (See also the discussion in Garcia et al. 2008). The ro- 
tation rate inferred below this radius is an a priori knowl- 
edge introduced via the regularization term. Then, 1 g mode 
-the candidate £=2 and n=-3- has been added with a typic 
error bar of 7.5 nHz (dotted line in Figure 7). In this case a 
constrain in the core is fixed by this mode and a better rate 
inside the core is obtained. When 8 g modes (£=l and 2 with 
n ranging form -2 to -6) are added (all with 7.5 nHz error 
bars) the agreement is even better (dashed line in Figure 7) 
than in the previous case with smaller error bars. The profile 
is still poorly determined and a higher number of modes or 
smaller error bars are needed to improve the information on 
the profile. Nevertheless, the resolution kernels in the latest 
case are much better determined as compared to the case of 
adding one single g mode. In this sense, for the deepest ker- 
nels, there are less wiggles and the secondary lobes become 
lower than with only 1 g mode (for more details see Mathur 
et al. 2008). For the deeper layers, below 0. 12i? , by adding 
more g modes, the secondary lobes are decreased by a factor 
of two and some of them even disappear. While the upper 
limit of the main lobe goes up to ~O.25i? in the first case, 
it decreases down to ~O.15i? in the second one. Thus, we 
can say that these resolution kernels are better localised. 
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Fig. 8 Correlation rates between the reference model (seismic) 
and another model with the same rotation profile. These models, 
described in Mathur et al. (2007), are: std93 (blue), std05 (red), 
tacho93 (light green), tacho05 (dark green), no diffusion (black) 
and model S (pink). The rotation profiles used have a differential 
rotation in the convective zone and a rigid rotation from 0.7 down 
to n (= r/Ro): 0.1 (crosses), 0.15 (asterisks) and 0.2 (diamonds). 
The rotation rate below n is constant: 3Q ra d (dotted line), 5Q ra a 
(dashed line) and 10^ ra d (dotted-dashed line), where £l ra d is the 
rotation rate in the rest of the radiative zone (433 nHz). Figure 
taken from Garcia, Mathur and Ballot (2008). 



4.2 Constraints for the solar structure 

The situation for the case of the structural physical param- 
eters is much better than for the solar rotation profile. In 
this case, for example, the profiles of sound speed and den- 
sity are rather well determined inside the solar interior even 
inside the solar core (e.g. Basu et al 1997; Couvidat, Turck- 
Chieze and Kosovichev 2003; Mathur et al. 2007) thanks 
to the information carried out by the £=0 modes that go 
through the deepest layers of the Sun. However, it would be 
better if we could add some constraints coming mainly from 
the core because the p modes -due to the increasing sound 
speed towards the interior- spend a much smaller fraction 
of time in these central regions than in the external layers. 

While we wait for the advent of the first individual g- 
mode frequencies, we can already check if recent measure- 
ments of the APi separation could be used to distinguish 
between models containing different physical processes in- 
side. Theoretical works have shown that the most recent 
solar models have APi differences of less than a minute 
(Mathur et al. 2007; Zaatri et al. 2007) which are too small 
to be directly detectable. However, instead of using the sep- 
arations themselves, we could utilize the correlation of the 
reconstructed waves (Garcia et al. 2007; Garcia, Mathur & 
Ballot 2008). By doing so, we have more sensitivity to the 
models and we could start to constrain some physical pro- 
cesses in the core of the Sun. 

The main uncertainty will come from the dynamics as 
the position of the m-components of each mode will de- 
pend on the rotational splitting. For this reason, to test the 
method, the correlations between a reference model (the 



saclay seismic model (Couvidat, Turck-Chieze and Koso- 
vichev 2003)) and different rotation rates (the same used 
in Garcia et al. (2007)) were calculated. Figure 8 shows 
the results of the comparison with 6 solar models includ- 
ing the model S (Christensen-Dalsgaard et al. 1996), stan- 
dard models (named std) and others that mainly differ on 
three physical inputs: the microscopic diffusion (all includ- 
ing this process but the one named no diffusion), the tur- 
bulent horizontal diffusion in the tachocline (named tacho), 
and the abundances of chemical elements. The models us- 
ing the old abundances (Grevesse and Noels 1993, named 
93) are clearly separated from those using newer ones (As- 
plund, Grevesse and Sauval 2005, named 05) and the ob- 
solete model without including the microscopic diffusion. 
This result is in agreement with previous analyses looking 
directly at the p-mode frequency differences computed for 
different models. However, in our analysis it is impossible 
to distinguish between models using very close physics like 
the standard models or those including the turbulence treat- 
ment in the tachocline. 

5 Conclusions 

In the present work we have shown that the longer GOLF 
time series confirm the results concerning the g-mode re- 
search. On one hand, there are still interesting features in 
the GOLF spectrum around 220.7 /iHz that are at the detec- 
tion limit of a quadruplet above the 90% confidence level. 
On the other hand, the amplitude of the pattern attributed 
to the signature of the dipole gravity modes has increased 
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and it contains more relative power than before. The highest 
peaks on the periodogram of the power spectrum could be 
interpreted as coming from APi, AP 2 and their harmonics. 

The study of the VIRGO data have revealed the pres- 
ence of one of the components of the quadruplet detected by 
GOLF in the same frequency bin. The analysis of the tempo- 
ral evolution of this signal in all the experiments of VIRGO 
has shown that this signal has been above noise level since 
the beginning of the mission and during the last 10 years. 
The SPM data has the highest signal-to noise ratio which 
is 3 orders of magnitude higher than the variations found 
in the same region for the house-keeping data, in particular, 
the temperature of the detector. 

The same analysis performed on velocity experiments 
has shown a similar behavior but with smaller signal-to- 
noise ratio. It is important to notice that the GONG data 
is a ground-based facility completely independent to SoHO 
and it would be very important if this preliminary result can 
be confirmed. 

Finally, we have illustrated how 2D inversions of the 
equatorial rotation profile can be improved by including the 
information given by 1 or several g modes and also how the 
use of the asymptotic properties of the g modes could be 
useful to put new constraints to the solar models. 
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